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We study conservative finite-mass corrections to the motion of a particle in a bound (eccentric) 
strong-field orbit around a Schwarzschild black hole. We assume the particle's mass /i is much smaller 
than the black hole mass M, and explore post-geodesic corrections of 0{^/M). Our analysis uses 
numerical data from a recently developed code that outputs the Lorenz-gauge gravitational self- 
force (GSF) acting on the particle along the eccentric geodesic. First, we calculate the 0(p./M) 
conservative correction to the periastron advance of the orbit, as a function of the (gauge-dependent) 
semilatus rectum and eccentricity. A gauge-invariant description of the GSF precession effect is 
made possible in the circular-orbit limit, where we express the correction to the periastron advance 
as a function of the invariant azimuthal frequency. We compare this relation with results from fully 
nonlinear numerical-relativistic simulations. In order to obtain a gauge- invariant measure of the GSF 
effect for fully eccentric orbits, we introduce a suitable generalization of Detweiler's circular-orbit 
"redshift" invariant. We compute the 0{n/M) conservative correction to this invariant, expressed 
as a function of the two invariant frequencies that parametrize the orbit. Our results are in good 
agreement with results from post-Newtonian calculations in the weak-field regime, as we shall report 
elsewhere. The results of our study can inform the development of analytical models for the dynamics 
of strongly gravitating binaries. They also provide an accurate benchmark for future numerical- 
relativistic simulations. 

I. INTRODUCTION AND SUMMARY 

We recently reported [ij a numerical code for computing the gravitational self-force (GSF) experienced by a mass 
particle set in motion around a Schwarzschild black hole. The code takes as input the two parameters of an eccen- 
tric geodesic {p and e, representing relativistic generalizations of the Keplerian semilatus rectum and eccentricity, 
respectively) and returns the various components of the GSF, evaluated in the Lorenz gauge, as functions along the 
geodesic orbit. Our current code does not take into account the back reaction from the GSF on the orbital motion; 
describing the radiative evolution of the orbit remains an important long-term goal of the self- force program 
strongly motivated within the context of experimental gravitational- wave physics |3|. However, GSF data from our 
current code (and from GSF codes by others d-Q, so far restricted to circular orbits) already give access to some 
interesting new physics. Specifically, they allow a quantitative description of some conservative post-geodesic effects 
associated with the finiteness of the particle's mass. This description is exact at linear order in the particle's mass fi 
(within a controlled numerical error, which in practice can be kept very small). 

Several such post-geodesic effects have already been analyzed. In Q Detweiler calculated the conservative 0(/i) 
correction to the value of the gauge-invariant "redshift" parameter {dtJdT, with t and r being the particle's coordinate 
and proper times, respectively; see below) along a circular orbit. In [7[ we derived the conservative O(^) shift in the 
frequency of the innermost stable circular orbit (ISCO). Most recently, in (with Damour) we computed the O(^) 
precession effect of the conservative GSF at the circular-orbit limit. These results were all successfully tested against 
the analytic predictions of post-Newtonian (FN) theory in the weak-field regime [1, [1, @, and, moreover, used to 
constrain some of the yet-unknown high-order FN parameters [sj-fllj . 

Indeed, GSF calculations already provide a valuable source of strong-field calibration data for various analytic 
models of the dynamics in binary systems. Huerta and Gair [l^l used circular-orbit GSF data from Ref. [l^ to assess 
the influence of conservative GSF corrections on the long-term phase evolution in astrophysical extreme-mass-ratio 
inspirals. Lousto et al. [ij, [l^ used the ISCO-shift result to inform an "empirical" formula (based also on results 
from fully nonlinear numerical simulations and on FN information) for the remnant masses and spins in binary black 
hole mergers. The value of the ISGO shift was also utilized by Favata [l^ as an exact reference point in an exhaustive 
study of the relative performances of various FN methods. And, in a recent collaboration with Damour Q we used 
GSF data for slightly eccentric orbits to constrain, for the first time, the strong-field shape of some of the analytic 
functions appearing in the effective one body (FOB) description of the conservative binary dynamics. In all these 
examples GSF data are valuable because they give a handle on an hitherto inaccessible, "remote" region of the binary 
parameter space, namely the domain of small mass-ratios and small separations. 

So far, analysis of the conservative effects of the GSF has been restricted to circular orbits or to orbits with an 
infinitesimally small eccentricity. Here we take a wide step forward, and study the conservative effects of the GSF 



2 



for orbits with large eccentricities ("large" in the sense of "not necessarily small" — in practice our code can handle 
eccentricities up to about 0.7; see below). In doing so we exploit the full capacity of our eccentric-orbit code for the 
first time. 

Eccentric geodesies of the Schwarzschild geometry are characterized by two frequencies: an azimuthal frequency ^l^p 
and a radial frequency fir, associated with the periodic motions in the Schwarzschild coordinates </? and r, respectively 
(the two frequencies will be defined more precisely below). Unlike in the analogous Keplerian problem where fl^p — ilr, 
here one always finds flip > fir, giving rise to what we interpret as periastron advance. The conservative piece of 
the GSF induces 0{^) shifts in the values of fl^p and fir and hence in the value of the periastron advance. Our first 
goal in this paper is to compute this GSF-induced correction to the periastron advance, denoted A(5, and we do so in 
Sec. [ni We express A(5 as a function of (the GSF-perturbed values of) p and e — a convenient parametrization of the 
orbit which is commonly adopted in perturbative studies (e.g., |17Hl9| V Unfortunately, although the shift A(5 itself 
is gauge invariant (in a sense that will be made precise in Sec. IIII A[) . the parameters p and e are not invariant, so 
our numerical results for A(5(p, e) must be interpreted with caution: they are valid specifically in the Lorenz gauge. 
A sample of our results for A(5(p, e) is presented in Tables ITTl and IIIII and plotted in Fig. [1] Notably, we find that the 
correction A(5(p, e) is negative for any p,e, i.e., the Lorenz-gauge GSF always acts to reduce the rate of periastron 
advance. 

The above calculation accounts only for the effect of the conservative piece of the GSF, and it assumes that 
its dissipative piece has been "turned off" . Of course, in the physical problem the frequencies fl^p and fir evolve 
dissipatively, and over one radial period they drift by an amount of 0(/i), which may well be comparable in magnitude 
to the conservative shift in these frequencies. One should therefore ask whether the conservative shift A5 has any 
physical relevance (notwithstanding the gauge issue mentioned above). We will suggest an interpretation of our results 
that makes them meaningful even in the presence of dissipation: The periastron advance in the physical, dissipatively 
evolving system may still be defined through the accumulated </3-phase between two successive periastron passages 
(minima of the orbital radius r); we will observe that this "true" periastron advance is given through 0(/i) by the 
conservative advance AJ calculated for a certain "average" geodesic. This observation makes it possible, in principle, 
to incorporate the 0(/i) precession effect in an approximate model of the orbital evolution (such as the one introduced 
recently by Gair et al. in Ref. 2]) based on our results for IS.5. 

It is clearly desirable to have at hand a gauge-invariant measure of the 0(/i) conservative effect — to allow, for 
example, a meaningful comparison with results from FN theory. What is required, more precisely, is a gauge- 
invariant relation, that expresses a nontrivial gauge-invariant quantity encoding the 0(/x) effect as a function of two 
other gauge-invariant quantities parametrizing the perturbed orbit. A natural gauge-invariant parametrization of 
the orbit is provided by the pair \ fl,„. r2r-}[38j. but it is not immediately obvious what one might choose for a third 
nontrivial invariant. In Sec. IIIII we will propose such a third invariant, which is constructed from the GSF (and the 
regularized metric perturbation) along the orbit, and encodes the 0(/i) conservative effect in a nontrivial way. Our 
gauge-invariant quantity, denoted (U)^, represents a natural generalization of Detweiler's "redshift" invariant: it is 
defined as the r- average of dt/dr over a radial period, which is also equal to the ratio between the radial period 
measured in time t and the radial period measured in time r. We then compute the conservative GSF correction 
A (U)^ as a function of the (GSF-corrected) frequencies fl^p and fir- Our numerical results for the gauge-invariant 
relation A {U)^ {flip, f^r) are displayed in Table ITVl Preliminary comparison of these results with FN expressions in the 
weak-field regime shows a very good agreement. A detailed comparison with FN theory will be presented elsewhere 

In the above description we have (for simplicity) left unmentioned a certain subtlety related to our choice of a time 
coordinate, and we come back to this point now. It is known [2l| that the monopole piece of the Lorenz-gauge 
metric perturbation has a peculiar behavior at r oo, which, however, can be cured using a simple "normalization" 
of the Schwarzschild time coordinate, in the form i — > t = (1 + a)t, where a is a certain 0{fi) constant which depends 
on the particle's orbit [see Eq. below]. The time t (unlike t) refiects the asymptotic fiatness of the perturbed 
spacetime, and at r ^> oo it coincides with the natural time coordinate used in FN studies. Since one of our main 
goals here is to provide useful data for FN comparisons, we choose to work with the time coordinate t throughout 
our analysis. In particular, we define the orbital frequencies with respect to time t, not t (and call them fl,p and fir 
to reflect this), and we similarly define (U)^ through dt/dr, not dt/dr. We note, however, that the quantity AS is 

insensitive to the time coordinate chosen (as it depends only on the ratio fl^p/flr = fl,p/flr). 

In Sec. lIVI we return to consider the periastron advance, this time focusing on the circular-orbit limit (e — ^ 0). This 
limit defines a one-parameter family of orbits, each of which nonetheless characterized by two invariant frequencies 
flip and fir (the latter associated with an infinitesimal e-perturbation of the circular orbit). This allows for a gauge- 
invariant description of the GSF precession effect in this case, e.g., by expressing 5 as a function flp. In previous 
work we have already analyzed the precession effect of the GSF for slightly eccentric orbits. Here we reproduce 
our results in a form that allows us a direct comparison with recent results from fully nonlinear numerical relativistic 
(NR) simulations by Mroue et al. ^] . We illustrate the feasibility and potential benefit of such comparisons in order 
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to motivate more detailed study. 

To help readers navigate through this work we give, in Table HI a key to our main symbols and essential notation. 
The table can serve as a quick reference guide for readers wishing to use our numerical results without delving into 
the technical details of their derivation. Throughout this work we use metric signature — I — h+ and "geometrized" 
units with c — G — 1 (with the mass of the central black hole providing a natural unit for both time and distance) . 
Physical units can be restored by multiplying all distances by G/c?, all frequencies by (? jG^ etc. 



symbol 


meaning 


where defined 


11 M 


particls's mELSS, blcick hole's mass 




n 
H 


jjj f A/f J small mass ratio 




t r 9 u> 


Sdiwarzschild/Lorenz-gauge coordinates 




peons 

^ OL 


conservative piece of the GSF (oc /i'^) 




Xn 


value of a quantity X at the geodesic limit (/i — ¥ 0) 




AX 


the Oi^fJi) conservative GSF correction to X 


e 0- Ens (l2(l (ISTIl 


p, e 


Qpmimi"iiQ "rppi"iiTn ( npT A/f\ ^nH pr'ppntvir'ii"^;" 

OClllilCLL no 1 Cl^ L Lllli I L/Cl IVl J CXil\-L Ctj^jCll Li. H^l L y 


Rn (TmIi 


X 


orbital radial phase ("true anomaly") 


Eq. (O 


E, L 


particle's energy and angular momentum per /i 


Eq. IHl) 




radial and azimuthal orbital frequencies 


Eqs. (dZl), dUl) 


T 


radial period (in time t) 




T 


radial period (in proper time r) 




6 


fractional periastron advance (dissipation ignored) 


Eq. (O 


^true 


fractional periastron advance (dissipation included) 


Eq. dlZl) 


i 


normalized ("asymptotically flat") time coordinate 


Eq. (O 


X 


a quantity X redefined with respect to time t 


e.g., Eq. dZO]) 




generalized redshift invariant 


Eq. (O 


To 


radius [r value) of a circular orbit 





TABLE I: Essential notation and key to main symbols. 



II. PRECESSION EFFECT OF THE GSF 



A. Preliminaries: geodesic motion 



Let us begin by recalling a few results from the theory of geodesic motion in Schwarzschild spacetime. Throughout 
our presentation we use a subscript '0' to distinguish "geodesic" quantities from their GSF-perturbed counterparts; 
thus, for example, Uq and u" will denote the tangent four-velocity associated with the geodesic and GSF-perturbed 
orbits, respectively. Timelike geodesies in Schwarzschild geometry constitute a 2-parameter family, parametrized by 
two constants of motion: the specific energy Eq = —Uto and specific angular momentum Lq = UipQ. Here Uao = gapUQ, 
where gap is the background Schwarzschild metric with mass M, and we have assumed (without loss of generality) 
that the geodesies lie in the "equatorial plane", 9 — Tr/2 (hence tig = 0). We will be interested in the subset of 
bound geodesies, i.e., ones confined to a radii range < r < for some Tq > AM ("periastron") and > 
("apastron") (there are no bound timelike geodesies with < AM). The pair {tq ^Tq ) can be used as an alternative 
parametrization of these geodesies, which we call eccentric. The special (one-parameter subset of) eccentric geodesies 
with = Vq are called circular. 

For eccentric geodesies it is convenient to define the (adimensionalized) semilatus rectum pq and eccentricity cq 
through 

PoM_ +_ PqM 

i + eo i — Co 

and reparametrize the orbits using the pair (po,eo). The structure of the {po,eo) parameter space is described, e.g., 
in Sec. ILA of [l| (see in particular Fig. 1 therein). Each eccentric geodesic has a unique value of (poiCq) within the 
range < Cq < 1 with pq > 6 + 2eo. Geodesies located along the separatrix pq = 6 + 2eo are marginally unstable (and 
will not remain bound under small perturbations of cq and/or po). Stable circular orbits have eo = and pq > 6. The 
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special geodesic with (po, cq) = (6, 0) (corresponding in the eo-po plane to the intersection of the separatrix with the 
eo = axis) is the ISCO. 

In terms of po and eo, the geodesic specific energy and angular momentum are given by 



En — 



PoiPo - 3 - e^) 
The orbital radius varies in time according to the formula 



(po-2-2eo)(po-2 + 2eo)]'/' ^ PqM 

' " ^ 7 O 2M /2 • ^ ) 



roix) = -r- , (3 

1 + eo cos X 



where the radial phase x ("true anomaly") is related to the coordinate time t — <o(x) through 



dto ^ Mpg / (po-2-2eo)(po-2 + 2eo) 

dx (po - 2 - 2eocosx)(l + eocosx)^ V po-6-2eoCOSx 



(4) 



with initial condition io(x = 0) = < where t corresponds to a periastron passage [i.e., rolxit )) = Tq]. The 
azimuthal phase of the orbit, ip — (paix), can be obtained as a function of x by solving 



dipo I PQ j-g^ 



dx V - 6 - 2eo cos x 

with initial condition ip{x = 0) = ip^ , where tp^ is the (/3-phase aX t — . We note that both io(x) a-nd <y5o(x) a^'S 
monotonically increasing functions. 

The f-period of the radial motion (i.e., the i-time interval between two successive periastron passages) and the 
frequency associated with it are given by 

To= -T^ dx, ^ro = 27r/ro, (6) 
Jo dx 

where the integrand can be expressed as a function of x using Eq. The azimuthal phase accumulated over one 
radial period To, i.e., $o = V'oix = 27r) - ipo{x = 0), is given by 

<^o - ^ = 4 ( l\ ) ellipK ( ) , (7) 

Jo dx \po~6-2eoJ \po-6-2eoJ 

where we have substituted from Eq. ([5]) and where ellipK(7) = /p^^^(l — 7 sin^ x)'~^^^dx is the complete elliptic integral 
of the first kind. The azimuthal frequency fl^po is defined as the t-average of {dip/dt)o over a radial period, which may 
also be expressed as 

rJ^o = *o/To^l}.o X [$o/(2^)]. (8) 

It can be shown that, for any given < eg < 1, $0 is a- monotonically decreasing function of po, with $0 ^ 27r as 
Po 00. Hence $0 > 27r for all eccentric orbits, with the excess $0 ~ 27r representing periastron advance. We shall 
denote the average periastron advance per radian by So- 

So = $0/(2^) - 1 > 0. (9) 

From Eq. ([8|) we then have the relation 

n^o/^rO=''i-+So- (10) 

It is instructive to examine some asymptotic characteristics of ^o- For a fixed eccentricity we have the large-po 
expansion 

Sipo » 1) - 'Spo' + ^(54 + 3e2)p-2 + Oip^^), (11) 
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and we note that cq only enters i^o at sub-leading order — and even there the effect of eccentricity is suppressed due to 
the smallness of the term Sbq compared with 54. We also observe that Sq does not vanish at the circular-orbit limit; 
rather, we find 

-5(60 « 1) = VPo/iPo - 6) - 1 + 0{el). (12) 

Note that at the ISCO (where flro = 0) Sq diverges as ^ {po — 6)^^^^. Lastly, consider the asymptotic behavior near 
the separatrix, po — 6 — 2eo <C 1 (for cq > 0): 

I f3 + eoY^^ , f 64eo 



In . (13) 

TT \ 2eo / Vpo - 6 - 2eo/ 

Note the logarithmic divergence at the separatrix limit. We see that the form of divergence of at the ISCO limit 
(POj eo) (6, 0) depends on the direction (in the Pq-cq plane) from which this limit is taken, illustrating the singular 
nature of the ISCO point. 

It should be noted that the pair of frequencies {Q.rO,^kpo) does not constitute a one-to-one parametrization of 
eccentric geodesic orbits in Schwarzschild geometry. The transformation (po,eo) (f^rO,^^<po) [and also {Eo,Lo) 
(f2rOi ^i/>o)] becomes singular along a certain curve in the parameter space, well outside the separatrix (see Fig. [3] in 
Appendix |^ ; one finds that for each orbit lying on the inner side of the singular curve there exists a "dual" orbit 
on the outer side, which is physically distinct but has precisely the same frequencies (rirOj^i^o)- This little-known 
property of the parameter space (we were unable to find any reference to it in the literature) is further discussed in 
Appendix 



B. The GSF-perturbed orbit 



Now consider a pointlike particle of mass /i <C M moving in the Schwarzschild geometry. At the limit /i — s> the 
particle's orbit is a geodesic of the background geometry, and we assume here this geodesic belongs to the class of 
eccentric geodesies discussed above. When the finiteness of /i is taken into account, we say that the particle experiences 
a GSF (cx /i^), and the motion is no longer strictly a geodesic of the background geometry. The GSF has a dissipative 
effect, which removes energy and angular momentum from the system (via gravitational radiation) and gives rise to 
a gradual inspiral. The GSF also has a conservative piece, which affects, for example, the precession rate of the orbit. 
The splitting of the GSF into dissipative and conservative pieces is formally defined in terms of the retarded and 
advanced metric perturbations (see, e.g., [1^); in Ref. il] we describe how each of the two pieces can be constructed 
from the full GSF in practice, taking advantage of the particular symmetries of Schwarzschild geodesies. We will 
proceed here by making the (nonphysical) assumption that the dissipative piece of the GSF has been "turned off" , 
and that the particle is moving under the influence of the conservative piece of the GSF alone, denoted F^°''^^. We 
shall come back to discuss the dissipative effect in the next subsection. 

We shall assume that the orbit remains strictly bound under the small perturbation caused by -F'c'ons (this is only 
allowed because we are ignoring dissipation), and denote the values of the new, "perturbed" radial turning points 
by r~ and r"*". (We note that the values of r^,r+ are gauge dependent, as are the values of many other perturbed 
quantities we define below. We will describe the construction of gauge-invariant quantities in the next section; in the 
meantime, for concreteness, we may assume that all perturbed quantities are given in a particular gauge — e.g., the 
Lorenz gauge.) We then define p and e via 

_ pM + pM 

T =1-—' 1 ' 

1 + e 1 — e 

and use the pair (p, e) to parametrize the perturbed orbit. We also define 

E = -ut, L = u^, (15) 

where Ua — gapu^ with being the 4-velocity tangent to the perturbed orbit, taken to be normalized with respect 
to the background metric: gapu'^u^ = —1. Note that E and L are not conserved along the orbit, and do not have the 
interpretation of energy and angular momentum. From symmetry, the orbit remains equatorial even under the effect 
of the GSF: = 0. 

We further assume that the perturbed orbit remains periodic in time, with radial period T interpreted as the 
perturbed value of Tq (that such a periodic solution to the perturbed equations of motion exists will be confirmed 
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below by explicit construction). We then choose a radial phase parameter x so defined that the radius along the 
perturbed orbit varies according to 

pM 

r = r[x) = — , (16) 

1 + e cos X 

in analogy with the geodesic formula ([3|) (we use the same symbol x s-s in the geodesic case for mere notational brevity; 
this should not lead to confusion). The effect of the GSF on the shape of the orbit is then encoded in the relations 
t = t{x) ^-iid = '^{x) along the perturbed orbit, which depend explicitly on the GSF (see below). The perturbed 
period and radial frequency are obtained through 

(it 

T= —dx, 2TT/T, (17) 

Jo dx 

and the perturbed accumulated phase $ and azimuthal frequency are given by 

^-fjpX, n, = ^/T. (18) 

Finally, the conservative GSF-perturbed value of the periastron advance reads 

5^<P/{2n)-l = n^/flr~l. (19) 

In what follows we will obtain an explicit expression for 5 as a function of p, e (which will, of course, involve F^°"-^). 
We will then define the 0(/i) "post-geodesic" conservative correction to 6 as 

AS = d{p,e)-Soip,e), (20) 

where So{p,e) is the "background" value obtained from Eq. ([9|) with Eq. ([7]), replacing (po,eo) — > {p,e). Notice 
that our definition of the ^^0{fi) correction" employs only the perturbed parameters (p, e). We are not considering 
here the alternative definition A6 = S{p,e) — (5o(po,eo), which, of course, differs from at the leading order, 
0(/x). The motivation for considering the correction A rather than A will become clear in the next section, where 
we reparametrize the orbit using two gauge-invariant quantities — the frequencies fi^ and fii^ — and calculate the A 
correction to a third invariant, expressed as a function of the two perturbed frequencies. This A correction will then 
be truly gauge invariant, unlike the corresponding A correction, which would depend on the background frequencies 
ilro and r2(^o- We remind, however, that, in our present discussion, neither A(5(p, e;po7 eo) ^lor A(5(p, e) constitute 
gauge-invariant relations. This is because the orbital parameters p, e themselves are gauge-ambiguous. In this section 
we will be calculating A(5 specifically in the Lorenz gauge. 

C. Conservative GSF correction to the periastron advance 

Recalling Eq. let us write 

^ = ¥'x(x;P,e) = <^xo(x;P>e) -I- Av3^(x;p,e), (21) 

where the "background" quantity '-Px^hiW^^') is obtained from Eq. (O replacing (pojCo) (p^e). From Eq. we 
then have 



2^ 7o 



¥'xo(x;P,e)dx - 1 



1 



1-n 



A^^{x;p,e)dx, (22) 



and, identifying the term in square brackets as Sq(p, e), we obtain from the definition in (j20p 



2 /■2Tr 



AS{p,e)^—l Aip^{x;p,e)dx. (23) 







To obtain an expression for Aip^{x^P^ g), let us start by writing 



^ ^ eL{x)\sinx\ , . 

'^^ rdx pM[E^{x)~V{r{x),L{xW/^' ^ ^ 
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where hereafter an overdot denotes d/dr, with r being proper time along the perturbed geodesic. Here we have used 
'■P = g'^'^u^ = r~'^L, substituting for r{x) from Eq. the factor dr/dx was evaluated using Eq. again; and 

to evaluate r we used the normalization of the perturbed 4-velocity, gapu^-u^ = —1, giving P ~ E"^ — ^(^(x); ^(x))j 
with the effective potential 



y(r,i)= 1-— 1 + -^ . (25) 



The "background" quantity Lpy^Q{x]p,e) in Eq. (|2T|) can be obtained from Eq. (f24|) by fixing p, e, x [hence also fixing 
r{x)] and replacing E{x;p,e) — ?> Eo(j>,e) and i(x;p, e) LQ(ji,e) [where i?o(p, e) and Lo{p,e) are the geodesic 
functional relations obtained from Eq. with (po, eo) (p, e)]- The GSF correction At/?^ — ip^ — ip^Q can therefore 
be obtained at 0(/i) by considering the A-variation of 93^ in Eq. p4| with respect to i? and L: 



p T dL 



AL. (26) 



Here the partial derivatives are taken with fixed p, e, x (hence also fixed r) and evaluated at (Eq{p, e), Lq{p, e)), and 
we have introduced AE{x;p,e) = E{x]p,e) — Eo{p,e) and AL{x;p,e) = L{x]p,e) — Lo{p,e). With Eqs. p3| and 
the task of calculating A(5 thus reduces to obtaining the 0(/i) corrections AE and AL. 
To obtain AE and AL we need to consider the t and tp components of the particle's equation of motion. These 
read, respectively, fiiit = — F™"** and fiii,^ — F™"**, or, with the definitions of Eq. (ITSI) . 

^lE = /iL = (27) 

ft is useful to think of F^°'^^ and F™"*' as (periodic) functions of x along the perturbed geodesic. We define the two 
0(/i) functions 



(28) 



rt-r dr 

S{x)^-f^-'J^ Fr^ix')—dx', C{x)^f^-'j^ Fr'ix')^dx'. 

The factor dr/dx i^i the integrands need only be evaluated at leading order in fj, [since the expressions in Eq. 
are already 0{n)], and we can therefore use for this purpose the geodesic expressions given earlier. Writing dr/dx ~ 
{dr/dt)o{dto/dx) = (1 - 2M/ro)E^'^{dtn/dx) and then substituting from Eqs. (0), Q and g]), we find, at leading 
order in /i. 



dr Mp^f^ 



P 



dx (1 + e cos x)^ y p — 6 — 2e cos x 



(29) 



Note we are using here the perturbed parameters (p, e) instead of {pojCq), which is allowed since interchanging 
(p, e) O {po,eo) affects the expressions in Eq. ([25]) only at 0{fi^). In Eq. (PSI) it is sufficient, at our working order, 
to evaluate the GSF components along the background geodesic (po,eQ). Given numerical data for the conservative 
GSF along the geodesic (the kind of data provided by our code [l|), the functions £{x) and C{x) can be evaluated 
numerically for any x through 0{ijl). In fact, it is sufficient to evaluate these two functions for < x < tt, given the 
"reflection" symmetry £{x) — ^^(27r — x) and C{x) — C{2-k — x) [see Ref. [H, where we explain that the conservative 
i, components of the GSF have the antisymmetry property Fi°'^^{x) = — Fj™"''(27r — x) and similarly for and 
note that the factors dr/dx in Eq. (pS)) are reflection-symmetric]. 
Integrating Eqs. (l27t . we now obtain 

E{x) = E{G)+£{x)=E^ + AE{Q)+£{x), 

L{x) = m + C{x) ^Lo + AL{0) + C{x), (30) 

where E{Q) and i(0) are the values of E and L at the periastron (x = 0), and the 0(/i) quantities Ai?(0) = 
E{0;p, e) — Eo{p, e) and AL(0) = L{0;p, e) — Lo{p, e) are the shifts in the values of E and L at the periastron due to 
the conservative GSF. We identify 

AEix) = AEiO)+£{x), 

AL{x) = Ai(0)+£(x). (31) 

The values of the shifts AE{Q) and AL(0) are found in the following manner. At the periastron (r = r~, x = 0) 
and apastron (r = r"*", x = tt), where f = 0, it follows from the normalization of the perturbed 4-velocity that 

E\Q) = F(r-, L(0)), E\^) = V{r+,L{7r)). (32) 
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The 0(/i) piece of these relations (holding p, e — hence also — fixed) is given by 

E^E{0) = (^l_^)^Ai(0), (33) 

EAEin) = (l - — j (^ryi^^W' (34) 

where AE{Tr) and AL{tt) are the 0(/i) shifts in the values of E and i at the apastron. These two are related to the 
shifts at the periastron through 

AE{7t) ^ AE{0)+£{tt), (35) 

AL{tt) = AL{0)+C{n). (36) 

Solving the four equations simultaneously for AE{Q), AL{0), AE{7t) and AL(7r), with the substitutions 

r~ — Mp/ (1 + e) and r+ = Mp/ (1 — e), we obtain 

^^^^^ ^ Mp-^s'- e^)^'^ - e)'(p - 2 + 2e)S£(^) - Sin)] , (37) 



AL(0) = 



1 



4e(p - 3 



[(l-e)2(p-2 + 2e)/:(7r)-£(7r)/i?] 



with 



B 



Eo{p, e)AfV Mpy^[{p - 2)2 - 4e2]i/2 ' 



(38) 



(39) 



where we have substituted for Eo{p,e) and Lo{p,e) from Eq. With these expressions, the functions AE{x) and 
AL{x) of Eq. (|31l) are now fully specified given the parameters p, e and the GSF functions £{x),C{x)- 

Finally, substituting for AE and AL from Eq. (1311) in Eq. (j26p . we obtain the rather unwieldy but explicit result 



p(p-3-e2)V2[(p_2)2-4e2]V2 
e2(p — 6 — 2ecosx)''/^ 

~Me2(p-6-2ecosx)^/2 



4cos2(x/2) 



sin^ X 



(1 



e)2(p- 



2 + 2e)Cin) [p{l + e^) - 2(1 + 3e2) + 2e(p - 3 - e2) cosx]/:(x) 



4cos2(x/2) 



sin^X 



(40) 



Recall Ai5 is obtained from Aip^ as prescribed in Eq. (|23l) . Notice in that equation that, since A6 and A(/3^ are of 0{p), 
we are allowed to replace (p, e) — >■ (poj eo) in their argument; this only results in an error of 0{iJ?) which we neglect in 
our treatment. We also note that Aip^^ has the "reflection" symmetry mentioned above, i.e., A(p^(x) = A(p^(27r — x)- 
This allows us to fold the integral piece j"^^ over onto . With these modifications, we write our final result for AS 
in the form 



A5(po,eo) 



1 



Av3x(x;Po,eo)dx, 



(41) 



where Aiy9p^(x;po, eo), given in Eq. (|40)) with the replacement (p, e) {po,eQ), only contains quantities which are 
evaluated along the background geodesic (po, gq). 

Let us summarize the construction of A6. One starts with (numerical) data sets for the conservative GSF compo- 
nents i^™"'^ and -F"™"^, evaluated as functions of x along the background geodesic {pq, cq) between x = a-nd X = 
These data are integrated via Eq. (|28| to yield the indefinite integrals £{x) and C{x), which in turn are used to 
construct the function A(p^{x) via Eq. (j40|) . A second integration [Eq. (|4T1) ] finally produces A6{po,eo)- Note in this 
procedure the GSF data are integrated over twice. It is possible to simplify this "integral of integral" structure by 
integrating by parts in Eq. (|41l) . We find in practice, however, that this does not lead to significant simplification in 
the numerical implementation, and we hence prefer to leave our working formula for A6 in its form (1411) . 

Inspecting Eq. (j40p . one may be worried about the apparent divergence of the terms oc sin~^ x at x = 0,7r and of 
the terms oc cos~2(t^/2) at x = tt. In fact, the function Aip^ is perfectly regular at these two points (and anywhere 
else within the domain < x ''')■ To see this, it suffices to notice the following. First, at the limit x ~^ have 
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E^C cx x^, since the two functions are even in % and they obviously vanish at x — ^ [from their definition in Eq. (|28l) . 
recahing F™"'^ and are regular functions of x along the orbit]. Hence f/sin^ x ^-^d Lj sin^ x have finite limits 

as X — 0. Second, at the limit x — we have the even Taylor expansion £ — £{n) + ^(x ^ TT)^£"{n) + 0{x — ti')'* 
[recalling the reflection symmetry £{x) — ^(^ ~ x)]: &n.d similarly for C; substituting these expansions in Eq. (|40|) 
reveals that the limit x ~^ of this expression is finite as well. Nonetheless, a direct numerical implementation of 
Eq. (PT|) does require some special care near x = 0, tt. In Appendix [B] we describe our method for dealing with this 
technical subtlety in practice. 

D. Numerical results 

Our numerical algorithm for calculating the GSF is described in detail in Ref. T| , where we also discuss the various 
sources of numerical error and our method for estimating the error bars on the numerical data. Our code returns 
the GSF as a function of x along a given geodesic defined by the input parameters po,eo. Our code can handle 
eccentricities in the range < eo ^ 0.7 and po values between the separatrix and ~ 100; outside this range our current 
algorithm is less effective and the computation burden becomes prohibitive (see [l| for details). We should point out 
that even in the above "workable" part of the parameter space, there are narrow "stripes" in the eo-po plane which 
are currently inaccessible to our code. These inaccessible stripes correspond to orbits for which the frequency ratio 
ilipo/rtro is close to a small integer (mainly 2 or 3). For such "resonant" orbits, our algorithm for calculating the 
dipole piece of the metric perturbation (which is based on a frequency domain method) becomes ineffective. Hence, 
for example, the orbits with {pQ, eo) — (8, 0.1) and (8, 0.2) have Ot^o/^^ro very close to 2, and are currently inaccessible 
to our code. We have discussed this problem in Ref. jlj] and are hoping to address it in future work. In the current 
study, however, we will content ourselves with simply avoiding the problematic orbits. 



po eo 


q-^A5 




Po 


eo 


q~^A5 


q~'AS/5o 


Po 


eo 


q-^AS 


q-^AS/So 


6.1 0.02 


-146(2) 


-20.7(2) 


7.5 


0.1 


-1.614(4) 


-1.298(3) 


10 


0.1 


-0.2392(3) 


-0.4111(5) 


6.2 0.05 


-57.0(2) 


-11.71(5) 


7.5 


0.2 


-1.788(3) 


-1.411(2) 


10 


0.2 


-0.2471(2) 


-0.4230(4) 


6.3 0.1 


-41.9(1) 


-10.23(3) 


7.5 


0.3 


-2.134(3) 


-1.629(2) 


10 


0.3 


-0.2606(2) 


-0.4433(3) 


6.4 0.1 


-19.71(5) 


-6.12(2) 


7.5 


0.4 


-2.778(3) 


-2.015(2) 


10 


0.4 


-0.2807(2) 


-0.4731(3) 


6.5 0.1 


-12.28(3) 


-4.51(1) 


7.5 


0.5 


-4.070(3) 


-2.733(2) 


10 


0.5 


-0.3084(3) 


-0.5138(4) 


6.5 0.2 


-31.15(6) 


-9.36(2) 


8.0 


0.3 


-1.140(1) 


-1.101(1) 


12 


0.1 


-0.10011(8) 


-0.2415(2) 


6.7 0.1 


-6.45(1) 


-3.010(6) 


8.0 


0.4 


-1.347(1) 


-1.264(1) 


12 


0.2 


-0.102104(9) 


-0.24580(2) 


6.7 0.2 


-9.45(1) 


-4.062(6) 


8.0 


0.5 


-1.682(2) 


-1.516(1) 


12 


0.3 


-0.105563(4) 


-0.25322(1) 


6.7 0.3 


-27.03(3) 


-9.30(1) 


8.5 


0.1 


-0.6156(5) 


-0.7276(6) 


12 


0.4 


-0.110620(2) 


-0.264025(5) 


7.0 0.1 


-3.372(9) 


-2.024(6) 


8.5 


0.2 


-0.6497(9) 


-0.762(1) 


12 


0.5 


-0.11752(3) 


-0.27867(6) 


7.0 0.2 


-4.089(6) 


-2.359(3) 


8.5 


0.3 


-0.7106(8) 


-0.822(1) 


15 


0.1 


-0.03960(1) 


-0.13604(5) 


7.0 0.3 


-5.909(7) 


-3.157(4) 


8.5 


0.4 


-0.8052(9) 


-0.914(1) 


15 


0.2 


-0.03996(4) 


-0.1371(1) 


7.0 0.4 


-11.99(1) 


-5.516(5) 


8.5 


0.5 


-0.9461(9) 


-1.046(1) 


15 


0.3 


-0.04057(4) 


-0.1389(2) 


7.0 0.45 


-24.64(2) 


-9.764(8) 


9.0 


0.1 


-0.4299(4) 


-0.5861(5) 


15 


0.4 


-0.04149(4) 


-0.1417(1) 


7.0 0.49 


-127.30(9) 


-37.25(3) 


9.0 


0.2 


-0.4493(6) 


-0.6089(8) 


15 


0.5 


-0.04278(3) 


-0.1455(1) 


7.0 0.499 


-1279.4(9) 


-268.3(2) 


9.0 


0.3 


-0.4832(5) 


-0.6483(7) 


20 


0.1 


-0.01353(1) 


-0.06932(8) 


7.0 0.4999 


-12735(7) 


-2075(1) 


9.0 


0.4 


-0.5346(6) 


-0.7070(7) 


20 


0.2 


-0.01348(3) 


-0.0690(2) 








9.0 


0.5 


-0.6082(6) 


-0.7894(8) 


20 


0.3 


-0.01338(3) 


-0.0684(1) 
















20 


0.4 


-0.01327(3) 


-0.0677(1) 
















20 


0.5 


-0.01314(3) 


-0.0669(1) 



TABLE II: Numerical results for the conservative (Lorenz-gauge) GSF correction to the periastron advance per radian, AS, 
for a sample of po, eo values, q = fi/M is the small mass ratio. The last column in each table shows the GSF correction as a 
fraction of the "background" advance So. Parenthetical figures indicate the estimated uncertainty in the last displayed decimals 
(due to numerical error); thus, e.g., —19.71(5) stands for —19.71 ± 0.05. Note that the effect of the Lorenz-gauge GSF is to 
reduce the periastron advance {AS < 0). Note also, in the data for po = 7.0, the manifest linear growth of AS with the inverse 
separatrix distance (po — 6 — 2eo)~^; this behavior is discussed in the text. 

Tables HIl IIIII and Figure [T] display a sample of numerical results for AS over the parameter space of p and e. The 
tables also compare A6 to the 'background' (geodesic) advance Sq, as derived from Eq. ^ [with Eq. ^]. We can 
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make the following observations. 

1. The conservative post-geodesic correction AS (in the Lorenz gauge) is negative for any po,eo; it causes recession 
of the periastron and acts to reduce the rate of periastron precession. 

2. At large po: exhibits a very weak dependence on cq. This is reminiscent of the behavior of the background 
advance Sq; recall Eq. (fTTj) . 

3. At large po, the magnitude of the GSF correction AS seems to fall off faster with po than that of the background 
advance So; the latter [recall Eq. (ITT]) ] faUs off as oc 1/po (at fixed cq). Indeed, Fig. [1] suggests the large-po 
behavior 

A5(po»l)cxp-4^ (42) 
where the proportionality constant (^ — 2000/^/Af) is independent of eg. 

4. Near the separatrix, the quantity AS diverges as oc (po — 6 — 2eo)^^. 



Po 


eo 


{e/q)AS 


ie/q)A5/So 


6.301 


0.15 


-3.58(1) 


-0.424(1) 


6.401 


0.20 


-3.204(6) 


-0.4265(8) 


6.451 


0.225 


-3.073(5) 


-0.4293(7) 


6.501 


0.250 


-2.967(4) 


-0.4328(7) 


6.551 


0.275 


-2.879(5) 


-0.4368(7) 


6.601 


0.30 


-2.808(4) 


-0.4414(7) 


6.651 


0.325 


-2.750(3) 


-0.4466(5) 


6.751 


0.375 


-2.661(3) 


-0.4580(5) 


6.801 


0.40 


-2.629(2) 


-0.4645(4) 


6.851 


0.425 


-2.603(2) 


-0.4712(4) 


6.901 


0.45 


-2.582(2) 


-0.4782(5) 


6.951 


0.475 


-2.567(2) 


-0.4858(4) 


7.001 


0.50 


-2.556(2) 


-0.4937(3) 


6.110 


0.05 


-5.98(3) 


-0.659(4) 


6.160 


0.075 


-4.91(2) 


-0.625(3) 


6.210 


0.10 


-4.29(1) 


-0.607(2) 


6.260 


0.125 


-3.89(1) 


-0.598(2) 



TABLE III: Additional numerical data for AS, exploring the near-separatrix regime. The structure of this table is similar to that 
of Table im but for convenience we have multiplied here the values of AS and AS /So by the separatrix distance e = po — 6 — 2eo. 
The entries in the upper part of the table correspond to e = 0.001, and those in the lower part to e = 0.01. 

The last of the above observations is particularly interesting. Recall that the background advance Sq diverges at 
the separatrix only logarithmically with po — 6 — 2eo [for fixed eo > 0; see Eq. ([131) ]. This means that, for any given 
mass ratio q, small as we wish, the GSF correction AS becomes comparable in magnitude to the background advance 
i5o at some separatrix distance po — 6 — 2eo 0(q), and would in fact become a dominant effect (giving rise to a net 
recession of the periastron!) at distances smaller still. Inspecting the data in Table [III we find, for example, that with 
q — 1/2075 the conservative GSF effect becomes comparable in magnitude (and opposite in sign) to the background 
effect for a near-separatrix geodesic with (po,eo) = (7,0.4999). In such a situation one would no longer trust our 
perturbative treatment, and the leading-order GSF would seem to be of a limited utility. 

It should be explained immediately, however, that we do not expect the above situation to realize itself in actual, 
inspiralling systems. In the physical system, radiation reaction "smears" the transition regime across an interval 
Ar = 0{q^/^) [24, 25], over which the orbit evolves quickly and is highly nongeodesic. For relevant (small) mass ratios 
q, the above condition po — 6 ~ 2eo ^ 0{q) ^ Ar sets the orbit to lie well within this transition regime, where the 
quantity AS, defined along closed geodesies, is no longer physically meaningful. We may ask, conversely, what the 
magnitude of AS/Sq is at the onset of the transition regime, i.e., for pq — 6 — 2eo ^ q"^^^. The data in Table HIl indicate 
that this magnitude may range from about 10~^ for q = 10^^ to about 10^"* for q = 10^^, and is therefore always 
quite small in the relevant range of mass ratios. 
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6 7 8 9 10 15 20 6 7 8 9 10 15 20 



FIG. 1: Plot of some of the numerical data presented in Table HIl The left and right panels show, respectively, the absolute and 
relative GSF corrections, AS and AS /So (divided by the mass ratio q = fi/M), as functions of po for a variety of eccentricities 
eo. The data are shown on a log-log scale, and we have in fact plotted —AS and —AS /So since AS itself is negative (while So is 
positive). Note AS depends very weakly on eo at large po, where it appears to fall off with a power-law oc Pq'^. The background 
advance So falls off as ~ 3p^^ [recall Eq. The dashed line in the left panel is a reference line AS = — 2000g/po- 



E. Effect of dissipation 



So far we have been making the (unphysical) assumption that the particle experiences no dissipative forces and 
that, consequently, the orbit remains strictly periodic. In reality, of course, the dynamics is nonconservative and 
the orbital parameters (e.g., E and L) drift secularly in time under the effect of the dissipative piece of the GSF. 
In a regime where the drift is sufficiently slow ( "adiabatic" ) we may still define the periastron advance by referring 
to the ip-phase accumulated over the period of time between two successive periastron passages (minima of r), 
notwithstanding the fact that two such periastra would occur at slightly different radii [if the radiation-reaction 
timescale is E/E ~ 0{M'^/fi), which we assume here, the two successive periastron radii would differ by a small 
amount cx /i]. Let us denote by Stmc the "true" periastron advance, defined between two periastron passages of the 
actual evolving orbit, taking into account all 0{fi) effects of the GSF, both conservative and dissipative (below we 
give a more precise definition of (5tiuc)- Here (and in Appendix [Cjl we will argue that, through 0{^), the quantity Jtruo 
is, in fact, given by the "conservative- only" advance S, calculated along a suitable "average" geodesic orbit. 

Let us make this statement more precise. First, we need a notion of "slowly evolving" p and e. For our purpose, it 
would suffice to concentrate on a particular radial cycle of the evolving orbit, between two apastron passages at times 
(say) ti and t2 {t2 > ti). We denote the radius of the (full GSF-perturbed) evolving orbit by f(i), so that fi = f{ti) 
and f2 = ^(^2) are two consecutive apastron radii of the true orbit. We assume that at times ti and t2 the true orbit 
is tangent to periodic orbits with parameters (pi, ei) and (p2, 62), respectively, which are geodesies perturbed by the 
conservative piece of the GSF only (in Lorenz gauge, for concreteness). Through 0{fi), then, the "slowly evolving" p 
and e are given, for ti < t < t2, hy the linear interpolation 

p(t)=Pi + 7^^(i-ii) + 0(M'), (43) 

12 — ti 

and similarly for e{t). Note that terms quadratic (and higher) in t — ti are also higher order in /i, and so are any 
conservative GSF corrections to the term linear in i — ti — both types of higher-order corrections can be neglected in 
our discussion. 

We further assume that the solution to the equations of motion for the true, evolving orbit can be obtained [through 
0(/z)] from the corresponding "conservative" solution (i.e., assuming only the conservative piece of the GSF is at play) 
via the simple replacement (p, e) — > [p{t),e(t)]. Hence, in particular, we assume [recalling Eq. (|T6l) ] 



1 + e(t(x))cosx 

where x is a certain parameter along the evolving orbit, and the relation t{x) is obtained (in principle) by replacing 
(p, e) — [p{t),e{t)] and x ~^ X in the "conservative" expression for dt/dx and integrating with the initial condition 
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t{xi) = ^1- Here xi is the value of x at the first apastron, determined from the conditions xi ~^ for /x — ;> 
together with {df / dx) j^^j^^ = through 0(/L^) (we determine the value xi explicitly in Appendix [C]) . Note the form 
f — r[x;p{t),e{t)], which reflects the two-timescale dependence of the orbital radius: f depends on a "fast" variable 
X (or t), as well as on the slowly varying parameters p and e. Writing f{x) as in Eq. (j44[) amounts to assuming that, 
at any fixed time t = to (ii < to < ^2), the relation r = f[x;p(tQ),e(tQ)] describes a "conservative" orbit of constant 
parameters p = p{to) and e = e(to), tangent to the true, evolving orbit. 
In a similar fashion, we write for the (^-phase along the evolving orbit 

^^^^ix)-f'Ax;pitm,mx))l (45) 

where the right-hand side is obtained from Eq. (|2T|) by replacing {p, e) — > [p{t), e{t)]. The total (p-phase accumulated 
over one radial period of the evolving orbit {— ^2 — ^i) is then given by 

*true= r ^x[x;p{tm,mx))]dx, (46) 

''Xl 

where, recall, xi = x(^i)j and we also denoted X2 = x(^2)- The "true" periastron advance (per radian) of the evolving 
orbit between ti and ^2 is given by 

^truc = (27r)-l$truc - 1- (47) 

With these definitions, our claim is that the true advance dtme, which accounts at 0(/i) for both conservative and 
dissipative effects of the GSF, can be computed through 



Strue = S{p,e)+0{^i^), (48) 



where 



P=^{Pi+P2), e=i(ei+e2). (49) 

The quantity 6{p, e) on the right-hand side of Eq. (|48)) is the "conservative-only" advance, calculated along the 
conservative GSF-perturbed orbit with "average" parameters (p, e). It is not difficult to convince oneself of the 
validity of Eq. ; a proof of this relation is presented in Appendix O Note that the "average" geodesic for use in 
Eq. (I48p can alternatively be defined through the average values of E and L, or of any other pair of slowly varying 
orbital parameters [because the secular drift in any such parameters over a radial period would be linear in t through 
0(yu)]. Note also that, through 0{ii), p, e are the momentary parameter values at the periastron of the evolving orbit. 

We point out that the relation (j48l) is gauge invariant in an obvious sense, since we are considering here the full 
[0{fjP) + 0{fi)] periastron advance without breaking it up into "background" and "GSF-correction" pieces as we have 
done before. One may, of course, also define the (gauge-dependent) GSF correction in Jtruc via A^truc = (^truc — ^o(p, e), 
which, by virtue of Eqs. ([20]) and (|48p . would relate to the conservative-only correction AS via 

Adtr-uc ^A6ip,e)+0{fi^). (50) 

Equations (pS)) and (150]) suggest a way in which our results for the conservative-only periastron advance 6 may be 
relevant to actual, evolving orbits: One only needs to reinterpret the quantity (5(p, e) as the advance (per radian) of 
an adiabatically evolving orbit, defined [via Eqs. p6)) . (|47)) ] between two periastra whose average parameter values 
are {p, e). This interpretation [which, we recall, is only valid through 0{p)] may form a basis for comparison between 
our GSF precession data and results from fully nonlinear numerical-rclativistic simulations of binary inspirals, once 
the latter are available for sufficiently small mass ratios. This interpretation also suggests a way in which our results 
for the conservative correction AS might be incorporated into one of the existing approximate (PN/pcrturbative) 
frameworks for computing the orbital evolution in extreme-mass-ratio systems [3.ll8l[l9]. which so far do not account 
for the conservative effects of the GSF. (Any such application of our results will, however, need to deal cautiously 
with the gauge dependence of A6.) 



III. GAUGE-INVARIANT EFFECT OF THE CONSERVATIVE GSF 



The quantity AS(p, e) considered above is gauge dependent, and as such it is of a limited utility. Our ambition in 
this section is to devise a gauge-invariant measure of the 0{fi) conservative GSF effect on eccentric geodesies. As 
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already mentioned, while S itself is gauge invariant (as it is constructed from the two fundamental frequencies of the 
orbit, which are gauge invariant), the 0{fj,) functional relation AS{p,e) is not invariant since the parameters p,e are 
gauge dependent. A gauge-invariant parametrization of the eccentric orbits is provided by the pair {i}^,ilr), but the 
relation S{rttp,ilr) is trivial [recall Eq. (IT^ ] and gives no information about the GSF. What is required is a third, 
independent gauge-invariant quantity (call it G) which depends in a nontrivial way on the two invariant frequencies. 
The relation AG{ft^,ftr) = G{iltp,ftr) ~ GQ{iltp,ftr) [where, recall, G(-) and Go(-) denote the functional relations 
with and without the GSF, respectively] would then provide a gauge-invariant measure of the GSF effect. 

The quest for such a function G is motivated by the wish to identify a common reference point for comparison 
between calculations of the GSF held in different gauges [l^ , and also to facilitate comparison with results from PN 
theory. We therefore remain mindful that our G would need to be readily accessible to both GSF and PN treatments. 
In what follows we propose such a gauge- invariant G and use our code to compute the GSF correction AG{n^, fir), 
ready for comparison with results from other methods when these become available. 



A. "Physically acceptable" gauge transformations 

First, we need to make precise the meaning of gauge invariance in our problem. An 0(/i) gauge transformation 

^ - ^" (51) 
changes the physical (retarded) metric perturbation /i^^ associated with particle by an amount 

S(hal3 = £,a-t3 + C/3;a, (52) 

where a semicolon denotes covariant differentiation with respect to the background (Schwarzschild) metric. In the 
Detweiler- Whiting interpretation [l^, the GSF is exerted by a certain smooth perturbation h^p derived from h^p. 
This so-called "R-ficld" (which is defined only in the local neighbourhood of the particle) transforms under ^" in the 
same way as hap (22,, ■ 

kh% = io.;p + ip.o.. (53) 

It is usually desirable to work in a gauge in which the metric perturbation {hap or h^p) correctly reflects the underlying 
symmetry of the physical system. In the case of eccentric geodesies, an essential requirement is for the metric 
perturbation to respect the periodicity of the orbit, i.e., to exhibit a bi-periodic spectrum with fundamental frequencies 
VLr and ^l^. In particular, when evaluated along the geodesic orbit (say, as a function of proper-time r), the field 
h^p need be periodic with a r-period T corresponding to the radial i-period T (this forbids, for example, a spurious 
secular growth of the perturbation over time). We write this periodicity condition as 

h%[x;{r)]=h%[x^^{T + r)] (54) 

(for any proper-time r along the orbit), where x^ — x'^{t) describes the geodesic orbit. In the special case of a circular 
orbit, the condition (|54p is replaced with the requirement of helical symmetry [26]. We require the R-field in any 
"physically acceptable" gauge to satisfy (IM)) . 

In discussing gauge invariance we wish to restrict attention to the class of transformations ^" that take one "physi- 
cally acceptable" perturbation to another. Such transformations need to preserve, in particular, the periodicity of the 
metric perturbation. This can be achieved by requiring that the generators ^" themselves are bi-periodic functions 
of t (with frequencies fi^ and 0.^), so that, in particular, 

n^(^(r)]=n2:^(r + r)] (55) 

for any r. In fact, to this class of transformations we may add t-translations of the form ^" oc tSf , which are 
themselves nonperiodic but do not interfere with the periodicity of the metric perturbation; we shall come back to 
this type of gauge displacements later in our discussion. For now, however, let us restrict attention to the class of 
gauge transformations satisfying the periodicity condition ([55]) . We shall refer to any quantity defined along the orbit 
as "gauge invariant" if it is invariant under any gauge transformation within this class. j39j 

Let us consider, for example, the frequencies 51^ and ^l^, already described above as "gauge invariant". It is readily 
seen that the t-period T is formally invariant under periodic gauge transformations: Changing integration variables 
in Eq. ^ gives 
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(taking r = at periastron without loss of generality) , which, under the transformation ((511) , is modified by an amount 

6iT=- 1^ ^dr^-em+ew^o, (57) 

following from the periodicity condition ([55)) [we hereafter use ^"(r) as a shorthand for ^"{x!^{t))]. Hence, the radial 
frequency fir — 2t:/T too is gauge invariant. Similarly, we have for the azimuthal frequency 

n, = irj*. (58) 



giving 



s^n^ = -^£^dr = -T-'iem em = o, (59) 

by virtue of the gauge invariance of T and the periodicity condition (1551) . This confirms that fl^p too is gauge invariant. 
Note that the formal invariance of fir and is a direct consequence of the periodicity condition ([55]); in general, the 
frequencies will not remain invariant under gauge transformations that fail to satisfy (|55|) . 



B. Generalized redshift invariant 

Detweiler (soj first pointed out that, in the case of a circular orbit, the quantity U = (i.e., the conservative 
GSF-perturbed t component of the particle's four- velocity) is invariant under gauge transformations that respect the 
helical symmetry of the black hole-particle configuration. A possible physical interpretation of U as an observable 
measure of gravitational redshift is discussed in Ref. l4|. The gauge- invariant relation {/(fi^) was later utilized for 
comparing between GSF results in different gauges [lal3l| and between GSF and PN results [1, [13, [HI ■ 

Here we propose a natural generalization of the redshift invariant to the case of eccentric orbits. Our proposed 
invariant is simply the r-average of U over a radial period, which is also, more simply, the ratio between the t and r 
periods: 

The quantity (U)^ is obviously gauge invariant, since both T and T are invariant. (We could as well choose T as our 
third invariant; we prefer (U)^ because it is dimensionless and because is reduces to the standard redshift invariant 
in the circular-orbit limit.) 

As in (26j . we make two small adjustments to our invariant, to better fit it for comparison with PN calculations. 
First, we replace the proper-time r (which is defined with respect to the background metric .g^^) with the proper-time 
f defined with respect to the perturbed metric ga^ + h^p- The two proper times are related, through 0{fi), via 

+ (61) 



where 

1 
2 



i/^ = -h^u'^uP, (62) 



with the R-field perturbation evaluated at the particle. (Here the four-velocity may be defined with respect 
to either r or f; the difference would affect only at sub-leading order in yit, which we neglect in our treatment.) 
The R-field combination can be constructed directly from the full (retarded) metric perturbation using a certain 
regularization procedure, which we describe in Appendix [D| The appendix also details the numerical computation of 
in practice, using our code. 

The second adjustment is a normalization of the time coordinate i, motivated as follows. As first mentioned in (2l| 
and discussed in length in [1, [^,[1^, the physical metric perturbation component hu, in the Lorenz-gauge, has the 
peculiarity that it does not vanish at r — > oo but rather approaches a constant nonzero value (which depends only on 
the orbital parameters) . This behavior is entirely attributed to the static piece of the mass monopole perturbation, 
and thus the asymptotic value of hu does not depend on the angular direction even for eccentric orbits. To remove this 
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gauge artifact and facilitate comparison with PN theory (where a more suitable "asymptotically flat" time coordinate 
is used), we introduce the normalized time coordinate 



where a = a(p, e) is given by 



i^{l + a)t, (63) 



-^hu{r ^ «)). (64) 



Then, through 0{fi), we have g^f + h^^ — > — 1 as r — > oo, as desired. In the circular-orbit case one finds [21| 
a(e = 0) = ^JL[ro{ro — 3Af)]^^/^, where r = To is the radius of the orbit. For eccentric orbits we do not have an 
analytic expression for a, but its numerical value can be extracted from our numerical solutions for the Lorenz-gauge 
perturbation, using Eq. (|64|. Note that the transformation t ^ t amounts to an O(^) gauge transformation ([5T|) with 
= — aifJ" = This transformation from the Lorenz-gauge time t to the "asymptotically flat" time t does not 
spoil the periodicity of the metric perturbation, but its generator ^" does not respect the periodicity condition (j55p . 
As a result, the transformation does not leave {U)^ invariant; in fact, one finds 5^ (U)^ — a (U)^. 
With the above two adjustments, we redefine our generalized redshift invariant as 

1 f'^ di _ f 



where 



dt (66) 



are [through 0(/i)] the radial periods measured in time t and proper time f, respectively. 



C. The geodesic limit of {U)- 

At the limit ^ Q the quantities ([/)^ and (f/)^ coincide, and are given by 

where Tq = To(po,eo) and To — ToipoiSo) are the geodesic radial periods measured in t and r, respectively. The 
period Tq is computed via Eqs. (jH]), which, for easy reference, we reproduce here in explicit form: 



The proper-time period To is obtained by integrating dr/dx given in Eq. (|29p : 

\-l/2 



Toip, e) = 2AT,3/2(, _ 3 _ T ip-e - 2ecosxy 

To (1-Kecosx)^ 



-dx- (69) 



D. Gauge-invariant parametrization of the orbit 

As already discussed, we are aiming to replace the p, e parametrization of the orbit with a gauge-invariant one, 
based on the (perturbed) fundamental frequencies flri^ip- More precisely, we wish to work with the "normalized" 
frequencies, defined with respect to time t. These are simply related to the original i-frequencies via 

nr = {i-a)nr, n^ = {i-a)n^, (70) 

valid through 0{fi). Rather than using the frequencies rir,Clip themselves, we find it convenient to introduce a new 
pair of gauge-invariant parameters, denoted p and e, which are obtained by inverting 



27r ^ $o(p,e) 

ro(p,e)' ro(p,e) 



= = (71) 
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where Tq and $o are the geodetic relations, given respectively in Eqs. ([68]) and (O (with the replacement p, e — )• p, e), 
and, recall, Clr and Jl^ are the GSF-perturbed frequencies. The quantities p and e are natural gauge-invariant notions 
of "semilatus rectum" and "eccentricity", in much the same way that the standard quantity = (Mfi^)"^/'^ (see, 
e.g., fs']) is a natural gauge-invariant notion of radius in the circular-orbit case. 

It is important to recall, however, that the relations f2r(p,e) and r2^(p,e) in Eq. ([71]) are not bijective, and thus 
cannot be inverted without a suitable restriction of the domain — we remind the reader of the discussion at the end of 
Sec. Ill Al and in Appendix [^ The inverse relations p(r2r, ^ip) and e(f2r, ^^p) are nonetheless well defined in each of the 
domains p > Ps(e) and p < Ps(e) (cf. Fig.[3|) separately. In what follows we will assume that the domain [p > Ps(e) or 
p < Ps(e)] has been prespecified and that p{rir,Vl^) and e{(lr,(lip) are the uniquely determined values corresponding 
to that domain. 

Even in a suitably restricted domain, it is not possible to invert Eq. (j7ip in explicit form to obtain p(f2rj ^ip) and 
e(r2,.,r2^). However, given the GSF in a particular gauge, we may express the invariants p,e through 0(/Lt) in terms 
of the gauge-dependent parameters e,p, using the linear variation formulas 



e = e 



dflr 

de 



AO, 



dp 



de 



AQ^ + Oifi^). 



(72) 
(73) 



Here the partial derivatives can be evaluated by inverting the transformation matrix d(rir,Cl^)/d{p,e), which itself 
can be computed numerically (for given p, e) based on Eq. ([7T|) with ([68)1 and (O (replacing p, e -> p, e). The subscript 
'0' indicates that these partial derivatives are to be evaluated at the geodesic limit. The 0(/i) quantities AClr and 
Atttp are the GSF corrections to the corresponding frequencies, defined as 



AClr{p, e) = tlr{p, e) - ilro{p, e), Af2^(p, e) = tl^{p, e) — riipo(p, e). 



(74) 



Unfortunately, the inversion formulas (|72p and (I73p become meaningless along the singular curve p = Ps(e), where 
the transformation matrix d{t!,r, riy)/9(p, e) is singular. In the following analysis we will therefore "keep away" from 
parameter-space points that lie directly on the singular curve. 

To calculate p(p, e) and e(p, e) in Eqs. ([72]) and ([73]) we need explicit expressions for Ailr and Ail^ in terms of p, e 
and the GSF. These are given, through 0{ii), by 



Aflr — ~^^rO 1 



An,, 



a 



f) 

$0 



AT 



(75) 
(76) 



where AT(j), e) = T(ji, e) 
themselves given by 



To(p, e) and A<i>(p, e) = <i>(p, e) - ^oiP,e) are the 0{fj.) GSF corrections to T and 



AT 



A$ = 



2 / Aip^{x;p,e)dx^2TTAS. 
Ja 



(77) 
(78) 



The quantity Aip-^, recall, is the GSF correction to ip-^ = dip/dx', it was given explicitly in Eq. (|40|) in terms of e, p 
and the GSF quantities £{x) and £(x)- The quantity At^ is the GSF correction to 



t. = — 



dt 



The evaluation of At^ is similar to that of Aip^ (see Sec. Ill C| . The result is 

Mj55/2(p-3 - 6^)1/2 \{p-2- 2e)(2 + e + ecosx) 



At^ 



6^(1 + e cosx)2(p — 6 — 2e cosx)3/2 
p(p-3-e2)i/2[(p_2)2_4e2]i/2 



+ 



e2(l -I- ecosx)^(p — 6 — 2ecosx)3/2 



(1 



8cos2(x/2) 
ecos xY 



p - 3 - + (1 + ecosx)' 



sin^ X 



■ 2 

sm X 



(1 - e)2(2 + e + ecosx) 
8cos2(x/2) 



(79) 

Six) 

(80) 
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As with A</3^ [Eq. (|40|) ]. here too the evaluation of the expression near the turning point x = 0,7r is subtle. The 
practical solution proposed in Appendix [B] applies here as well. 

Let us summarize. Given a parameter-space point p, e and the GSF corresponding to that point, Eqs. (|72[) and (|73|) 
are used to construct the two quantities p(p, e) and e(p, e). The pair (p, e) constitutes a gauge- invariant parametrization 
of the eccentric orbits [on each side of the singular curve p = Ps(e)]. 



E. Conservative GSF correction to (U)- 

We are now in position to write down a gauge-invariant expression for the post-geodesic correction to {U}-: 

A{U)^^{U)4p,e)-{U)^„{p,e). (81) 



Here {U)^^ (PiS) represents the geodesic functional relation given in Eq. ([67]); it is to be calculated using Eqs. ((68)) 
and ([M)) with the arguments p,e replaced with p,e. Since both the function (C/)~ and the parameter pair p, e are 

gauge invariant, the 0{fj,) quantity A(C/)- provides a genuinely invariant description of the GSF effect. 

Our next goal is to express Eq. (j81|) in a workable form, i.e., as a function of the parameters p,e and the GSF. 
To this end, we expand each of the two terms on the right-hand side of Eq. ([81]) about its geodesic value {U)^q [p, e) 
through O(^). Starting with (C/)-(p,e) = r(p, e)/T(p, e), we obtain 

/ AT AT \ 

(C/),(p, e) (p, e) 1 + —— ~ —1— + 0{n'). (82) 

\ To(p,e) To{p,e) J 

Here AT and AT are the GSF corrections to T and T, respectively, which we shall give explicitly below in terms of 
p, e and the GSF. The second term on the right-hand side of ((8T|) is expanded in the form 



{U)roiP,^) = {U),,{p,e) + 



d{U)rO 



dp 



de 



(e - e) + 0{fi^), (83) 



where the subscript '0' indicates that the partial derivatives are to be evaluated at (p,e) — )■ {p,e). Substituting Eqs. 
dHU and dMl) in Eq. (|HT|) and using Eqs. (|721) and ([731) we then obtain, neglecting 0(^2) terms, 

A{U), - (C/),o (p, e) (^^^^ - 7^^)) - CriP, e)Anr - C^{p, e)A^^ (84) 

with 

The coefficients Cr and C^p can be computed (numerically), for any given p, e, from the appropriate geodesic expres- 
sions: The partial derivatives of (f/)^g are obtained using Eq. (|67l) with (jBH]) and (|69l) . and the partial derivatives 
d{p,e)ld(ytr,^i^ are computed as explained below Eq. ([75|. Note that in the expressions for Cr and we have 

allowed ourselves to replace (il^j^^') {yt^^Vl^) and (p,e) — > (p, e), which does not affect A(J7)- through 0(/x°). 

Finally, to be able to use Eq. (|84|) . we need expressions for the 0(/i) quantities AT, AT, Af^r and Afii^. The latter 
two have already been given above, in Eqs. (|75|) and ([76l) . As for AT, it follows from Eq. (|66l) that 



AT = aTo(p,e)-f AT, (87) 
with AT given in Eq. ([TT]) [with ^^\. Lastly, we have, recalling Eq. 

At = 2 / (Ar^ - iJ^r^o) ^X, (88) 
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where we introduced the notation 



_ dr 

with T-f^oip,^) being the geodetic hmit of r^, and At^(p, e) = T^{p,e) — T^o{p,e) 
given exphcitly in Eq. (|29p . and for the perturbed quantity we write 

dr/dx epAf| sinx|(l + ecosx)"^ 



^ dr/dr [E^x)-V{r{x).L{xWn 
proceeding using the method of Sec. Ill Cl to obtain 

Mp2[(p-2)2-4e2]i/2(p_3_e2) r 



(89) 



The background quantity t^q is 



(90) 



Ar-v 



6^(1 + e cosxYiv — 6 — 2e cosx)-^/^ 



4cos2(x/2) 



g(x) ' 

sin^X 



pi/2(p-3-e2) 
e^d) — 6 — 2ecosx)'^^^ 



(l-e)2(p-2 + 2e) /:(^) 
(l + ecosx)2 4cos2(x/2) 



— (p — 2 — 2e cos x) 



C{x) ' 

sir? X. 



+ 



Mp2[(p _ 2)2 - 4e2]i/2£:(^) _ pi/2(i _ e)2(p _ 2 + 2e)£(7r) 



(91) 



4e(l + ecosx)^(p — 6 — 2ecosx)^/^ 
I can again be used to assist in evaluating the last expression near the turning points 



The method of Appendix | 
X = 0,7r. 

Let us summarize the above construction. Our main result is expressed in Eq. (|84p . giving A(t/)^ in terms of the 
parameters p, e and the GSF. The various elements of Eq. ((5^ are constructed using Eqs. (I571) - (p5|) . ((75|) . ([75)1 and 
(|85|) - (|88|) . Note that since A(?7)^ is already O(^), in Eq. (|84| we may replace (p, e) — ^ (po,eo) and calculate the 
various ingredients of this equation based on GSF data evaluated along geodesic orbits. 

The value of A([/)- in Eq. (j84p does not depend on the gauge in which the GSF is given, and as such it provides a 
useful reference for comparison between results from different calculation schemes (PN or perturbative) . 



F. Numerical results 

Using Eq. ([M)) we have computed A(f/)- for a sample of e,p values; the results are displayed in Table HVl We 
have tested these results against weak-field analytic expressions through IPN, and found a very good agreement — this 
comparison will be presented in a forthcoming paper [23|. In the future, the data could also provide a basis for 
comparison with other calculations of the GSF (in whatever gauge) when these are available. 

In selecting the dataset for display in Table IIVI we have deliberately avoided parameter-space points which are 
located on or near the singular curve p — Ps{e), where the quantity A{U)- becomes singular. This indeterminacy is 
an unavoidable price to pay for using the gauge- invariant parametrization (fir, il^ip), which is ill-suited along Ps{e). 

It is worth clarifying a potentially confusing point regrading the interpretation of the data in Table IIVI The table 
labels the orbits by their (p, e)-values and not by their invariant frequencies (Clr, (l^) or the invariant parameters (p, e) 
associated with them. This does not mean that we are measuring the GSF effect with respect to the background 
quantity {U) ^^{p^e): we remind that the GSF correction A(C/)- is defined in an invariant way with respect to 
(C/)^o(P'^) — recall Eq. (I8ip . The p, e parametrization is adopted a posteriori for convenience, and it by no means 
compromises the gauge invariance of A([/)-. 

Finally, we should comment on how the information in Table ITVl could be used for comparisons with PN expressions, 
the latter being usually given in terms of the two invariant frequencies Vlr,CLip (or simple combinations thereof; see, 
e.g., [13) and not in terms of p,e. Suppose one has a PN expression for A{U)~, given as a function of Clr,Clip. 
One starts by extracting from the PN expression all terms through 0(/Lt) (holding Clr, fixed). One then subtracts 
the "background" [0(^°)] quantity {U)fQ{0,r,(l^), where, of course, {lr,Cl,p are the full, perturbed frequencies. The 
remaining [0(/i)] terms are explicit functions of VlrT^ip, which may now be replaced with the background values 
ilro, ^ipo [with the error from this replacement being of only 0(/i2)]. The resulting expression is readily evaluated for 
any given values of p, e from Table llVl and the result may be compared with the corresponding GSF value of A(C/)~ 
given in the table. This procedure has been used to facilitate the comparison to be presented in Ref. j20j . 
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IV. THE CIRCULAR-ORBIT LIMIT AND A TENTATIVE COMPARISON WITH NUMERICAL 

RELATIVITY 



The limit e ^> (with fixed p) defines a circular orbit of radius pM = ro{= const), which, however, still has two 
distinct frequencies associated with it: the azimuthal frequency $7^ (or Clip) is simply dip/dt (or d(p/dt) along the 
limiting circular orbit, and the radial frequency fir (or Cl^) is that associated with a slightly eccentric orbit resulting 
from an infinitesimal e-perturbation of the circular orbit. Hence, rather conveniently for us here, the circular-orbit 
limit defines a one-parameter family of orbits that are nonetheless characterized by two gauge-invariant frequencies. 
We may utilize the relation between these two frequencies (or between any two independent combinations thereof) as 
a simple gauge-invariant function for GSF studies. In Ref. we have already studied the circular-orbit limit, and 
extracted the 0(/i) gauge-invariant information embedded in the relation between the two frequencies (then used it 
to inform a comparison with PN-calibrated BOB models). 

In a recent numerical relativistic (NR) study Mroue et al. report measurements of the periastron advance in fully 
nonlinear simulations of slightly eccentric black hole inspirals. Data are provided for a sample of binary mass ratios in 
the range 1:1-1:6, with eccentricities of order a few xlO^^ (see Fig. 8 of [l^; how the eccentricity is defined in these 
simulations is discussed in Sec. II therein). It is interesting to examine the NR data against the predictions of our GSF 
calculation in the circular-orbit limit. Of course, we must exercise great caution in attempting such a comparison, 
and we note here three obvious caveats. First, even the "small" mass ratio of 1:6 is well outside the natural domain 
of validity of our 0{fj,) GSF approximation. Second, the NR simulation automatically accounts for dissipative effects, 
which are ignored in our GSF calculation. Our argument (Sec. Ill Ep that dissipation has a negligible effect on the 
periastron advance loses its validity when the mass ratio is not sufficiently small. Third, Ref. [22| makes no special 
effort to obtain very accurate precession data, since its main motivation lies somewhere else (it attempts to develop 
a method for eliminating spurious eccentricity in quasicircular inspirals). As a result, the statistical variance (due to 
numerical error) in the precession data of [23| is large, making a meaningful comparison difficult. The numerical error 
is particularly large for smaller mass ratios, which are, alas, most useful to us here. 

Nonetheless, we would like here to take advantage of this opportunity to make a first contact between the GSF 
and NR programs, if only to point to the potential of a mutually beneficial synergy between the two programs, and 
in order to motive further study. 

Mroue et al. plot the frequency ratio Cl^p/Clr versus the adimensionalized azimuthal frequency {M + iJ,)Cl^. Note we 
must interpret the NR frequencies as our "hat" frequencies (those defined with respect to the "asymptotically flat" 
time t) rather than our "Lorenz-gauge" frequencies Qtp, It is also crucial to notice that Mroue et al. adimensionalize 
the frequency using the total mass M + n and not (as customary in GSF analysis) the large mass M. A GSF expression 
for Cl^/Clr is readily obtained, through 0(/x), using Eqs. (3), (6), and (14) of Ref . Q: 



Clr \ To 



-3/2 r 



^ _ 6A/ ro(ro -3M) ^, 



ro{ro-3M) 
2^M 



(ro-3M)3/^ ■ 

1 /iMl/2r2 ^_ 



(92) 



Here is the radius of the limiting circular orbit, and the 0{fi'^) quantities , and F^ arise from the formal 
e-expansion of the conservative GSF components along the slightly eccentric orbit. Specifically, F^ is the r component 
of the conservative GSF along the strictly circular orbit, and Ff,-?"^ are associated with the 0(e) variation of the 
GSF. We refer the reader to i§\ for a precise definition of these GSF quantities. How these quantities are extracted 
in practice from the numerical GSF data is described in Ref. [ij. The azimuthal frequency is given through 0(/x) by 




1 - an — 



rg(ro-3M) 
2nM{ro - 2M) 



f: 



(93) 



[see Eqs. (3) and (6) of g], where Uo = a{e — 0) = fi [ro(ro - 3A/)] 



-1/2 



Hence, through 0{ij,), 



iM + fi)Cl^ = (A//ro)=^/2 



/i - Q!o — 



rliro-3M) 
2nM{ro-2M)' 



(94) 



In Appendix lEl we give numerical values for the GSF coefficients F^{ro), Fl{ro) and F^{ro). Using these values in 

Eqs. and we can obtain a numerical (parametric) relation between Cl^/Clr{— fJ^/ilr) and (M -I- ^)r2^ for 
any given value of /j,. We plot this relation in Fig. [2] for a sample of mass ratios q = fs/M . Superposed on these GSF 
"predictions" (which, of course, represent gross extrapolations beyond the natural mass-ratio domain of the GSF) we 
display the NR data points from Mroue et al. [l^]. Each NR data set, for a given g, comes from a single inspiral 
simulation; how the frequencies are extracted from the numerical data is explained in Ref. [23 |. 



20 




FIG. 2: (color online) Tentative comparison of GSF and NR data for the periastron advance of slightly eccentric orbits. Shown 
is the frequency ratio f^^^/fir as a function of Cl,p (adimensionalized using the total mass M + fi), for a variety of mass ratios 
q — fJ^/M between 1:2 and 1:6. The dashed line corresponds to a test particle {q = 0). Single data points describe results 
from NR simulations, reproduced here from Fig. 8 of Mroue et al. Solid lines are interpolated 0{q) GSF "predictions", 

calculated using Eqs. (|92p and (|94p with the numerical values of the GSF coefficients tabulated in AppendixlE] The horizontal 
scale of this plot roughly coincides with that of Fig. 8 of [22| for easy reference. Despite the manifest low accuracy of the NR 
data, this preliminary comparison is already rather instructive (as described in the text), and motivates further study. 



The NR data clearly resolve a nonzero precession effect, but they are not accurate enough to allow a detailed 
quantitative comparison with the GSF predictions. Nevertheless, one can make several tentative observations, (i) 
The NR and GSF data are in agreement on the "sign" of the post-geodesic precession effect: it is opposite that of the 
geodesic precession, i.e., the GSF acts to reduce the rate of periastron advance, (ii) The NR data for q = 1/2 is least 
noisy and perhaps most accurate. If we are to trust these data, we find that the GSF prediction "overestimates" the 
post-geodesic effect by about a factor 2. This is not unreasonable for a mass ratio as large as 1:2. We see how our 
comparison already starts to tell us about the sign and magnitude of uncalculated higher-order GSF contributions. 
Tentatively, it would seem that the 2nd-order GSF precession effect is opposite in sense to that of the Ist-order effect, 
(iii) A most meaningful comparison (relatively speaking) would have been provided by the 1:6 data, if not for the 
very large scatter of the NR points in this case. The 1:6 NR data seem roughly evenly distributed about the GSF 
curve, but it is not possible to make more definite statements. 

The above tentative comparison illustrates the potential benefits from synergic GSF/NR studies. Foremostly, it 
provides a strong two-way test of the results, because GSF and NR computations use highly independent methods. 
From the point of view of NR practitioners, the GSF predictions provide an accurate benchmark against which to 
assess the quality of the numerical simulations. From the GSF point of view, comparison with NR simulations gives 
access to valuable information about the effect of currently inaccessible high-order GSF corrections. We envisage 
a synergy between GSF and NR methods as a fast-track avenue to the modelling of the two-body dynamics in 
intermediate mass-ratio inspirals, which is currently beyond the reach of either method. 

Motivated by the above, we have recently initiated a collaborative study to pursue and exploit a more detailed 
comparison of the GSF and NR precession data. Preliminary new NR data by Mroue et al. are dramatically more 
accurate, and show a remarkable agreement with the GSF predictions at q — 1/8. We hope to report results from 
this study in a forthcoming paper [33| . 
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Appendix A: Isofrequency geodesic orbits 

We think there is a common belief that bound geodesic orbits in Schwarzschild (or Kerr) spacetime can be labelled 
uniquely by their frequencies. This belief is unfounded, and turns out to be false. To the best of our knowledge, this 
issue was (rather surprisingly) never addressed in the literature, so we briefly discuss it here. 

In Fig. Owe plot the level lines n;^o=const and r2rO=const over a portion of the po, eo parameter space of bound 
geodesies in Schwarzschild. Recall the region of parameter space with po > 6 + 2eo (and < eg < 1) corresponds 
to bound geodesies. The transformation {po,eo) — >■ (JlrOj^^ipo) becomes singular along the separatrix po = 6 + 2eo. 
More surprisingly, it also turns out to be singular along a certain curve po = Ps(eo), well outside the separatrix, which 
divides the parameter space into two disjoint domains, po < Ps{eo) and po > Ps{eo) (see Fig. [3|). At each point along 
Ps(eo), the i7^o=const level line is tangent to an r2ro=const level line, and the Jacobian matrix of the transformation 
(potSq) — s- (rirOj^^ipo) becomes singular. Note that the transformation {Eq,Lq) (pQ,eo) is perfectly regular across 
the entire parameter space of bound geodesies, so one cannot dismiss the above behavior as a mere peculiarity of the 
Po,eo parametrization. 




FIG. 3: (color online) Singularity of tlie transformation (po,eo) — > (f2^o,^^ro) for bound (eccentric) geodesies in Schwarzschild 
spacetime. The plot displays the level lines r2^o=const (dotted, red) and firo^const (solid, blue) over a portion of the po,eo 
parameter space. The diagonal po — 6 + 2eo (straight black line) is the separatrix, where QrO = 0; stable geodesic orbits exist in 
the region po > 6 + 2eo. In the domain shown, Q^o decreases monotonically with po while ilro increases monotonically with po 
(the sign of dQ,ro/dpo reverses further out to the right of the region shown). The thick (black) curve, po = p3(eo), is the locus 
of points where the fi^o and SlrO level lines are tangent to one another; along this curve the Jacobian of the transformation 
(po,eo) — > (Q.^o,Q.ro) vanishes. Each orbit left of the singular curve has a "dual" isofrequency orbit associated with it, which 
lies to the right of the singular curve. One such pair is indicated in the plot (black crosses), at (po,eo) = (6.3,0.05) and 
(6.59274, 0.27569). 

It can be easily verified that the presence of the singular curve po — Ps{so) results in the following: for each orbit 
belonging to the "left" domain [po < ps{eo)], there exists a "dual" (or "isofrequency") orbit in the "right" domain [po > 
Ps(eo)], which is physically distinct but has the same frequencies n^Q,^rO- For example, the orbit with parameters 
{PO: eo) = (6.3, 0.05) has the same frequencies as an orbit with parameters {po, eo) = (6.59274 . . . , 0.27569 . . .). This, 
of course, means that the frequency pair (fl^oT^ro) does not label the geodesic orbits uniquely. However, one may 
still use the (fi^Oj^^ro) parametrization separately in each of the "left" and "right" domains (as we do in Sec. IIIip . 
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It should be noted that all orbits on the "left" domain are deep within the zoom-whirl regime. The "least" zoom- 
whirling orbit in the left domain [the one corresponding to the intersection of Psiso) with eg = 0] has a frequency 
ratio J7,po/f^rO — 4.1, i.e., it executes over 3 whirls per radial period. It is also interesting to ask whether the range 
of dual orbits in the "right" domain extends to the weak-field regime (which would suggest fascinating astrophysical 
implications). It is clear, however, that this is not the case: the azimuthal frequencies of orbits in the left domain 
(and hence those of their "right" duals too) are confined to the range 

6.38-3/2 < ^.^^^^ ^ 4-3/2_ 

Readers familiar with the literature on radiation reaction in black hole spacetimes may find the curve Ps(eo) in Fig. 
12] here reminiscent of the "critical curve" shown in Fig. 3 of Cutler et al. [17,], which describes the locus of points in the 
Po, Co plane at which the radiative evolution of the eccentricity changes its sense (so that de/dp turns from negative to 
positive). Inspection of the two plots reveals, however, that the two curves lie sufficiently far apart in the parameter 
space to assume that there is no direct relation between them. It is still interesting to ask about possible anomalies 
in the behavior of a radiatively inspiraling object as it crosses the curve Ps{eo)- If the radiative evolution happens to 
drive the inspiral across Ps{eo) in a direction (in the po,eo plane) tangent, or nearly tangent, to the frequency level 
lines, we might expect a "hang up " ep isode during which the evolution of the two frequencies halts, or slows down. 
The information in Fig. 3 of Rcf. 17] suggests, however, that the radiative evolution through Ps{eo) proceeds in a 
direction roughly orthogonal to the frequency level lines. Any "hang up" effect should therefore be minimal. 

One may also ponder the possibility of finding resonant interaction effects acting between pairs of dual orbits around 
astrophysical black holes (e.g., between clumps of accreting matter). This intriguing possibility deserves exploration. 
A preliminary step would be to understand the dual behavior across the full (3-dimensional) parameter space of 
generic orbits in Kerr geometry. We have confirmed that a similar duality exists at least in the subspace of eccentric 
equatorial geodesies in Kerr [341] . 

Finally, we point to the fact that, in the context of the GSF problem, the relation between the two frequencies 
along the (GSF-perturbed) singular curve, fl^ = fl^psii^r), is a gauge-invariant one. As such, we can envisage it being 
utilized as a reference for comparison between different calculations of the GSF, and for a strong-field calibration of 
approximate analytic methods — in much the same way that the relation fl^{ilr) has been utilized in the circular-orbit 
limit . We hope to explore this possibility in future work. 



Appendix B: Treatment of numerical x-integrals near periastron and apastron 



The numerical implementation of formulas (pij) (for AS), ((77)) (for AT) and (|88l) (for AT) is somewhat subtle, due 
to the formal divergence of individual terms in the respective integrands A.ip^ , At^ and Ar^ at the two radial turning 
points, X = 0, TT. As we explained in the text, each of the full integrands is in fact smooth for all x, including at 
X = 0, TT. However, in practice, the divergence of individual terms require a special treatment at the turning points, 
which we describe here. In what follows we refer specifically to A6 (and Aip^^) for concreteness; AT and AT are 
treated in a similar manner. 

Consider first the periapsis, x = 0, which is more easily dealt with. The function Aip^{x) [Eq. (HO]) ] has the form 



A^x = /i(x) + /2(x)^ + /3(x) 
sm X 



• 2 ' 

sm X 



(Bl) 



where /i, /2, /a are certain functions of x which are regular (smooth) at x = {At^ and At^ have similar forms, with 
different /„'s). The functions £{x) and C{x) have even Taylor expansions at x = 0, with f (0) = £(0) = 0, and so the 
expression in Eq. (|B1|) is in fact perfectly regular at x = 0. Still, the factors £/ sin^ x a-^d C/ sin^ x pose a practical 
problem, because the numerical integration routine we apply to evaluate Eqs. (j4ip requires as input (also) the value 
A(y9^(0). We deal with this simply by writing 



f(x) 

hm — 2 — 

x-^o sin X 

hm^ 

x-^o sin X 



1 
2 

1 £_C 

2 dx^ 



1 

2^ 





dT\ 




V dx 


dx) 


x=o 



1 



x=o 





dT\ 




V dx 


dx) 


x=o 



(B2) 



where we have recalled the definitions of £ and C in Eqs. ((28)) . The factor dr/dx is given explicitly in Eq. ((29)) . and 
in the above equalities we have made use of the fact that it is an even function of x, so the term oc d'^r /dx^ vanishes 
at X = 0. The GSF derivatives in Eqs. (jB2p are readily evaluated numerically at x = from our GSF data sets. Then 
the numerical value of A(^^(0) is obtained by substituting from Eqs. (IB2p in the form ((Bip . 
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Now turn to consider the apastron, x = ti": where the situation is shghtly more involved. For this discussion, we 
note that A(p^ can also be written in the form 



ffi(x)fW -52(x)^(x) 



. 2 

sm X 



sin^ X 



where the g„'s are certain functions of x which are smooth near x = ti", where they satisfy 

51 (X) - 52(x) + 0{x - ^)', g3(x) = 34(x) + 0(x - Tr)^. 



(B3) 



(B4) 



(Ai^ can be written in the same form, with different functions having the same properties. The same also applies 
to At;^-, modulo an additive function of x which is however smooth at x = tt and hence of no concern to us here.) Since 
£ and C admit regular Taylor expansions about x = tt, of the form £ = £{tt) + 0{x ~ tt)^ and L = C{tt) + 0{x — tt)^. 



it follows that the expression in Eq. (|B3P is perfectly regular at x = and, in particular, the limit x ~^ of this 
expressions is finite. Here, however, it is not sufficient to obtain A(/3^(7r) via a Taylor expansion as we did for x = 0. 



An added practical difficulty is that, for values of x near tt, a delicate cancellation occurs in Eq. (IB3p between the gi 
and g2 terms and also between the 53 and g^ terms. This can give rise to large numerical errors due to roundoff if 
one attempts to evaluate Eq. (jB3P directly. 
To circumvent this problem, we introduce 



TTicons 



TT'Cons 



^£in)-£ix), 

^W-Ax), 



(B5) 



in terms of which Eq. (|B3|) becomes 



^Vx^hiix)£{TT)+g2ix) 



g(x) 

sin^ X 



h2{x)C{7r) + g^ix) 



SVC? X 



(B6) 



where hi{x) = [.9i(x) ~52(x)]/ sin^ X and /i2(x) = [33(x) ~ff4(x)]/ sin^ x are regular (smooth) at x = tt. This form no 



longer involves a delicate cancellation between different terms near x 
difficulty. To obtain A(^^(7r) we simply use 



TT and is thus free from the above numerical 



limiW 

x^-ir sin X 

hm^ 

x-i-T sin X 
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2 dx^ 
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2 dx^ 







dT\ 






K dx 


dx) 


X=Tr 
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2fl 





dT\ 


V dx 


dx) 



(B7) 



where the GSF derivatives are evaluated from the numerical data. 

Finally, to carry out the integral in Eq. (j4ip . we split the integration domain as 



A^=i 

TT 



7r/2 ^ 

A(p^dx + - 



tt/2 



A'fxdx, 



(B8) 



where we use the form (jBip [with (|B2p ] for the first integral, and the form (|B6p [with (|B7p ] for the second. The 
integrals for AT [Eq. ([77)) ] and for AT [Eq. ([88]) ] are evaluated in the same way. 



Appendix C: Relation between "true" and "conservative-only" periastron advance 

In this appendix we establish the relation (j48]), which states that the periastron advance <5truG of the physical, 
evolving orbit [as it is defined in Eq. (H71) ] is equal through 0(/i) to the conservative advance 6 associated with a 
certain "average" conservative orbit. 

We begin by reproducing Eqs. (|46|) and (j47|) here for easy reference: 

<5truc = (27r)-i$truc - 1, $truc = / p(i (x)) , e(t (x) )]dx- (CI) 

J XI 
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Here, recall, the integration is over a complete radial cycle of the slowly evolving orbit, from an apastron at x = Xi 
(t = ti) to the next apastron at x = X2 (t = ^2)- The values of xi and X2 are fomid from the conditions xi — tt 
and X2 — J' TT for — 0, together with r'(xi) — ^'(X2) = 0, where throughout this appendix a prime denotes d/dx- 
Using Eq. (|33]), these conditions give 



Xi = -TT + (^xbi^ei) + O(At^), 

X2 = TT + Sx{P2,e2)+0{fl^), 



with 



6xip, e) 



- e) +pe' 



pe 



(C2) 



(C3) 



Here it is sufficient to evaluate p' and e' at leading order [0(/i)], which may be done with the help of Eqs. (|43p and 
([4]). Recall (pi,ei) and (^2,62) are the orbital parameters at times ti and t2, respectively (assumed given). Since 6x 
is already 0(/i), in Eq. (jC2[) we are allowed to replace pi,p2 — > P and ei, 62 — ^ e, where p, e are the average parameter 
values defined in Eq. p9)) . p, e are also, through 0{^), the parameter values at the periastron. We thus have 



Xi = -TT + Sxip,e) + 0(p'^), 
X2 = Tr + Sx{p,e) + 0(p'^). 



(C4) 



Now let us formally expand p and e in x about the periastron through 0{ij), noting that the periastron value of x 
is 0(/i): 



P^P + p'{p, e)x + 0{^i'^), e = e + e'{p, e)x + 



(C5) 



where p',e' [each of 0(/i)] are evaluated at the periastron. We use this to expand the integrand in Eq. (jCip in the 
form 



X2 

dx 

'XI 

$1 + $2 



Vxix\P^ e) + 



dVxix.P. e) 



i9p 



de 



e'x + 0(/x2) 



(C6) 



where $1, $2 and $3 denote the corresponding contributions to $truc from the first, second and third terms in the 
integrand [neglecting the terms of 0(^^)]. 

Let us consider $2 first. Using Eq. (|C4|) we can express it as 



$2 = 



dp 



dp 



p'ttSx 



d^^{x;p,e) 



dp 



P'X + 0{i? 



(C7) 



The first two terms on the right-hand side may be evaluated at the geodesic limit, since they are multiplied by 
5x 0{fi). We note, recalling Eq. ([5]), that iy9^(x;p, e) is an even function of X: and so is its partial derivative with 
respect to p. Therefore, the first two terms on the right-hand side in Eq. (jC7p cancel each other through 0(/i). In 
addition, we note that the integrand in the third term is odd in x and it follows that the integral vanishes. Hence, 
the entire contribution $2 is 0(/i^), and similarly for $3: 



$2,$3 = 0(/i2) 

Concentrate then on $1. Using Eq. (jC4l) again we have 



^ [ip^{Tr;p,e) - (fix{-'^;P,e)]Sx+ / (px{x;P,e)dx + Oip"^), 



(C8) 



(C9) 



in which the first two terms cancel each other through 0(/z) by virtue of the aforementioned even parity of (p-^{x',P, e), 
and the integral is simply the quantity ^{p, e) associated with the conservative orbit with parameters p, e. Hence, 



$1 = $(p,e)-F0(/x2). 
Combining Eqs. jCl]), and (|C10|) we conclude 

$truc = Hp, e) + 0{n^), 
and substituting $truo = 27r((5truc + 1) and e) — 2tt[S{p, e) -I- 1] leads directly to Eq. 



(CIO) 
(Cll) 
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Appendix D: Mode-sum regularization of H 



In this appendix we describe the calculation of the R-field combination 



(evaluated along the orbit) that goes into the expression for our gauge invariant A{U)f [Eq. (|84p with ([55]) in Sec. lIII| . 
Our numerical code returns the tensor-harmonic modes of the retarded metric perturbation (not those of the i?-field 
h^^), and the construction of the function follows through a certain regularization procedure, which resembles 
the standard mode-sum regularization of the GSF. Here we derive the necessary "mode-sum regularization formula" 
for the quantity H^. Our analysis follows very closely the method of Refs. (SSmTf . although the details are much 
simpler in our present case. 
We start by recalling 

^a/3=/»a"-/^f/3, (D2) 

where ft,^^' is the full (retarded) metric perturbation (denoted simply /i^^ elsewhere in this work) and ft,^^ is Detweiler- 
Whiting's S-field [13 • Both h^^ and h^p are uniquely defined (in terms of a particular Green function) in a small 
neighborhood of the particle's worldline, with the field h^^ being a smooth {C°°) homogeneous solution of the 
linearized field equations. At a field point a; = in the vicinity of a point z = on the worldline, the S-field admits 
the local expansion [3^1 

= [I + 0{6x')] , (D3) 

where an overbar denotes trace reversal, Vci{x',z) is the four- velocity vector parallel-propagated from z to a; along a 
short geodesic section connecting the two points, e{x; z) is the spatial geodesic distance from x to the worldline, and 
Sx = X — z. The spatial distance e itself has the local expansion 

e = eo (^1 + ^) + OiSx% (D4) 

in which ep = {gap + UaUp)5x°'5x^ (with the metric function and four velocity evaluated at z) and Si is a cubic 
polynomial in 5x [given explicitly in Eq. (A5) of [35]]. Defining now the field H^{x) = ^h^pv^'v^, we find 

hH^) = I - ^^^.S^^So/?) v^vf^ = ^ + OiSx), (D5) 

which, using Eq. (jD4[) . gives 

H^{x) = ii - ^ + 0{Sx). (D6) 
£0 2eo 



Let us now introduce the field i?*""(x) = ^h^afi'^°'^^ ■ Then, noting v°'{x — ^ z) = u", we obtain using Eqs. (|Dip 
and ini]) 

= lim [i7f""(a;) - H^ix)] . (D7) 

As in the standard mode-sum regularization prescription, we formally expand both H^^^^{x) and H^{x) in spherical 
harmonics (on a 2-sphere of constant r, and write Eq. (|D7I) in the form 

oo 

= lim ^ [Hl^\x) - Hf{x)\ , (D8) 

where represent the total /-mode contribution (summed over m) to i/f^'i/S'^ j^g shown, based on the 

form of the local expansion in Eq. (ID6p , that each of the /-modes Hi {x) is continuous (though not differentiable) at 
cc — > z; and that, at large /, H^{z) admits an expansion of the form 

Hi'{z) = BH + CH/l + 0{l-^), (D9) 
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where the coefficients Bh and Ch depend on z but not on I. Since iJ^"''(a;) — H^{x) is a smooth function, the 
multipole sum in Eq. (|D8p must converge uniformly and faster than any power-law in and it follows that Hf^'^^{z), 
too, must admit the large-l expansion 

= Bh + Ch/1 + 0{l-^), (DIO) 



with the same Bj{ and Ch- This allows us to reexpress Eq. (jD8|) as a sum of two convergent series, using the familiar 
mode-sum form 

oo 

H^ = J2 -^Bh- Ch/1] - Dh, (Dll) 

1=0 

with 

oo 

DH = Y.[Hf^^)-^H-CH/l\- (D12) 

(=0 

The parameters Bh, Ch and Dh are akin to the standard self- force "regularization parameters" — which we have 
attempted to reflect in our notation. The absence of a regularization term cx / is related to the fact that H^, unlike 
the self force, does not involve derivatives of the metric perturbation. 

The derivation of the parameters Bh, Ch and Dh follows closely the method detailed in Sees. VII and VIII of 
Ref. (ssj . It is based on expanding H^{x) [given in Eq. (jD6|) ] in spherical harmonics, and then evaluating the Z-mode 
contribution at the limit x — ?• z. First, one notes that the terms of 0{5x) and higher in Eq. (ID6p can be discarded 
as they cannot affect the value of in Eq. (jD7[) . [The individual Z-modes of the 0{5x) terms may well be nonzero, 
but their summed contribution must vanish at the limit a; — > z.] We can therefore write 

= lim [ijf"'Ha;) - iJ^^^PP™'^(a;)l , (D13) 

where 

jjS,approx(^) = ii _ (D14) 

£0 Sep 

and consider the ^-mode contributions from ijS.a,ppmx jj^g^^gad of those from . Considering first the term oc 5*1, 
one readily shows based on a simple symmetry consideration (cf. Sec. VII. B of [3^1) that it yields no contribution to 
Hi '^PP'°^(^z). The validity of this consideration does not depend on the explicit form of 5*1 (but only on the fact that 
it is cubic in 6x), and it is the same consideration that leads one to conclude the vanishing of the 'C parameter in 
the self-force case. 

Hence, the sole contribution to Hf''^^^'^°^{z) comes from the term /i/eo in Eq. (|D14p . This contribution is easily 
evaluated in explicit form using the method of Sec. VII. C of |35|, and one finds, crucially, that it is /-independent. One 
therefore identifies this contribution with the parameter Bh [recall Eq. (ID9p ]. and moreover concludes that Ch — 0. 
Furthermore, since the above implies '^pp''™^^) _ Bh ~ Ch/1 ~ for all I, it follows from the definition in Eq. 
(|DT2l) that Dh = 0. Exphcitly, one finds, in summary, 

where, recall, ellipK(-) is the complete elliptic integral of the first kind. Equation (jDllI) thus reduces to 

oo 

H^ = Y,[Hr{z)-BH]. (Die) 

(=0 

To calculate along the orbit we begin by recording the values of the (Lorenz-gauge) tensor harmonic /-modes 
' generated by our time-domain code (the /-modes are continuous at the orbit, so these values are well defined). 
We then construct the full /-mode contributions iJ[""(z) (as functions along the orbit) using 

We use these modes as input in the mode-sum formula (ID16I) . which yields H^. In practice we compute numerically 
only the first 20 or so /-modes. We estimate the contribution from the uncomputed large-/ tail by fitting a power-law 
model {a/P + b/l^) to the numerical data, then add the estimated tail contribution to the mode-sum. In this we 
follow the same procedure as for the GSF, and we refer the reader to Ref. [l| for details. 
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Appendix E: Numerical values for the circular-limit GSF coefiicients 

Table IVl gives numerical values for the GSF coefficients F^' , F[ and appearing in Eqs. (|92)) and (|94|. for a sample 
of circular-orbit radii To . How these values are extracted in practice from the numerical GSF data is explained in Ref. 

We give these values here for the benefit of readers who wish to reproduce the GSF curves shown in Fig. [51 or 
wish to obtain similar curves for other values of the mass ratio q. 
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An important caveat is discussed in the last paragraph of Sec. Ill Al and in Appendix 1X1 

The technical notion of gauge invariance in the case of a quantity that is only defined along the orbit is subtly different 
from that of a field like h^^. The operator in Eq. H53p . for instance, measures the difference in the functional form of 
h^lj{x) as a result of the transformation (|5ip . i.e., the difference between the 'new' and 'old' fields evaluated at the same 
coordinate value (and hence, generally, at different physical points). On the other hand, for quantities like u" or Q^p, which 
are only defined along the orbit, one considers the difference between the 'new' and 'old' quantities evaluated at the same 
physical point (same value of r), which in general takes different coordinate values before and after the transformation. 
Ref. [U] discusses this point in detail. 
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TABLE IV: GSF correction to the generalized redshift invariant {U)-: sample numerical results. Each row in the table 
corresponds to a particular geodesic orbit with parameters (p, e) and associated frequencies ($7^0, f2<po), as specified. The 5th 
column gives the numerical value of the coefHcient a{p, e) [see Eq. (|64[l ] needed for converting the Lorenz-gauge time coordinate 
t to the "asymptotically flat" time t. The sixth column displays the numerical values of the gauge-invariant quantity A{U}-, 
computed via Eq. (f84|) . Parenthetical figures indicate the estimated uncertainty in the last displayed decimals. 
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q'^F[ X 10^ 




80 


0.29897884(1) 


0.58707883(1) 


-0.318667(3) 


50 


0.74494870(4) 1.4474199(1) 


-0.391467(4) 


40 


1.1428834(1) 


2.2055600(3) 


-0.429392(1) 


30 


1.96981742(2) 3.7606729(4) 


-0.4805562(2) 
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2.76576978(2) 5.2382103(4) 


-0.5136676(2) 
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-0.629220(3) 
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10.046259(3) 


18.50998(2) 


-0.643288(3) 
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11.551775(4) 


21.33090(3) 


-0.659417(3) 


10 


13.389514(2) 


24.89511(5) 


-0.679283(4) 
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15.637164(4) 


29.54408(9) 


-0.706774(4) 


8.5 


16.936718(5) 


32.4544(1) 


-0.725949(5) 


8.0 


18.357927(4) 


35.9162(2) 


-0.751553(2) 


7.5 


19.89068(1) 


40.1309(2) 


-0.78753(2) 
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21.49921(1) 


45.4273(2) 


-0.84103(2) 


6.5 


23.09361(1) 


52.3795(3) 


-0.92545(3) 


6.0 


24.4665(1) 


62.095(1) 


-1.0665(8) 



TABLE V: Numerical values for the circular- limit GSF coefficients appearing in Eqs. (|92p and ((94}. Vo is the radius of the 
limiting circular orbit. Parenthetical figures show the estimated uncertainty in the last displayed decimals. Data from this 
table were used to generate the GSF curves in Fig. [2l 



